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Abstract. We present a detailed analytical study of a paradigmatic scale-free network model, the Static 
Model. Analytical expressions for its main properties are derived by using the hidden variables formalism. 
We map the model into a canonic hidden variables one, and solve the latter. The good agreement between 
our predictions and extensive simulations of the original model suggests that the mapping is exact in the 
infinite network size limit. One of the most remarkable findings of this study is the presence of relevant 
disassortative correlations, which are induced by the physical condition of absence of self and multiple 
connections. 
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1 Introduction 

In the last few years a considerable amount of research ef- 
fort has been devoted to the study of a large array of natu- 
ral and man-made systems that can be described in terms 
of networks. In fact, systems as diverse as the Internet .1 , 
the World-Wide Web [2], collaborations networks 0(3; 
the web of sexual contacts |S] , foodwebs |H]|7j , protein in- 
teractions networks |HI , metabolic networks [Hj , and many 
others can be represented, at a certain level of approxi- 
mation, as networks or graphs ^H], in which vertices rep- 
resent the elementary units composing the system, while 
the edges stand for the relations or interactions present 
between pairs of elements. These kind of systems were 
in the previous century the subject of study of classical 
graph theory. Recently, however, the availability of large 
data sets and more powerful computer resources, together 
with the application of new statistical tools, has led to the 
development of a modern theory of complex networks 
IT^ , which is nowadays one of the most active fields in the 
statistical physics of complex systems. 

Despite of the wide variety of the systems considered 
in this field, some common characteristics seems to be 
present in almost all complex networks. Among them, the 
most remarkable is probably the fact that many real- world 
networks exhibit a fat-tailed degree distribution P{k). That 
is, the probability that a randomly chosen vertex has a 
number of emerging edges equal to an integer k (the de- 
gree of the vertex) has the form for large k 

P{k)^k-\ (1) 

where the degree exponent 7 commonly ranges in the in- 
terval 7 G [2,3] This suggest the presence of a het- 
erogeneous hierarchy of vertices, lacking a characteristic 
degree value, which has result in the common denomina- 
tion of scale-free networks fT^. Moreover, the presence of a 



scale-free degree distribution implies that the degree fluc- 
tuations are unbounded in the infinite network size limit, 
i.e. (/c^) 00 when N 00, which has a considerable im- 
pact on the behavior of dynamical processes taking place 
on top of the network. For instance, it has been shown 
that scale-free networks are extremely resilient to random 
damage |14lll5lll6j . while at the same time they are very 
weak in front to the spread of epidemic processes |17II18| . 

In addition to the degree distribution, it has been po- 
inted out that real networks are further characterized by 
the presence of degree correlations, which translate in the 
fact that the degrees at the end points of any given edge 
are not usually independent. This kind of degree corre- 
lations can be quantitatively expressed in terms of the 
conditional probability P{k'\k) that a vertex of degree k 
is connected to a vertex of degree k' |19II2(J) . From a nu- 
merical point of view, the presence of correlations can be 
conveniently studied by means of the average degree of 
the nearest neighbors of the vertices of degree fc, formally 
defined as ^Hl 

knn{k)=Y.^'P{k'\k). (2) 
k' 

A first classification of networks has been proposed ac- 
cording to the nature of their correlations ^T) . Thus, when 
knn{k) is a growing function of k, the network is said to 
exhibit assortative mixing, while a decreasing fc„,i(fc) func- 
tion is typical of disassortative mixing. 

The appealing evidence for the existence of scale-free 
networks has prompted the development of numerous mod- 
els, aimed at understanding the origin of fat-tailed degree 
distributions, or even the nature of degree correlations 
P^ . Scale- free network models can be roughly divided in 
two main classes: Growing network models, capitalizing in 
the original Barabasi- Albert model ^j, focus their ap- 
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proach on the evolution of the network, rather than on its 
structure. The key ingredient of these models consist in 
considering the network as a result of a growth process, 
in which new vertices and edges are sequentially added to 
the system following a prescribed set of dynamical rules 
(usually inspired in the preferential attachment or rich- 
get-richer paradigm J13'). On the other hand, static net- 
work models consider networks with a constant number of 
vertices N , among which edges are drawn following differ- 
ent probabilistic rules. In this sense, Ref. 23 proposed for 
the first time a static network model yielding a scale-free 
degree distribution, while Ref. ,20; (see also Ref. 24 ) pro- 
posed a general class of static network models, the hidden 
variables network models, which allow to develop a sys- 
tematic analytical formalism for this class of systems. 

One of the models belonging to the class of static net- 
work models that has recently attracted some attention, 
is the so-called static model (SM), recently proposed in 
Ref. The interest raised by the SM has a twofold 

origin. Firstly, its definition is very simple, and allows to 
generate large networks with any desired degree exponent 
with a reasonable amount of computer effort. Secondly, its 
use as a benchmark to check the properties of both scale- 
free networks and dynamical processes running on top of 
networks has become quite widespread lately 25,26,27, 
I28ll29j . In spite of this wide interest and use, however, lit- 
tle is known from an analytical point of view about the 
properties of the networks generated by the SM (apart 
from its scale-free nature), especially in what refers the 
nature of its possible degree correlations. 

In this paper we present an analytical treatment of 
the SM based in a mapping to a hidden variables network 
model. Using the hidden variables formalism, we are able 
to provide analytic expressions for the main properties of 
the SM, in particular for the degree distribution and the 
degree correlations, as measured by the fc„„(/c) function, 
showing a very good agreement with extensive numerical 
simulations of the original SM model. One of the most re- 
markable findings that we report is the presence of strong 
disassortative degree correlations in the SM for values of 
the degree exponents close to 2, in agreement with the the- 
oretical arguments put forward in Ref. pITj . The presence 
of this correlations indicate that the results of dynami- 
cal processes running on top of networks generated with 
the SM should be interpreted with great care, in order 
to discern the effects due to the scale-free nature of the 
networks from those related with the presence of intrinsic 
degree correlations. 

The paper is organized as follows. In Section |3 we re- 
view the definition of the SM model, as well as some of 
the its properties, that can be derived by using simple 
qualitative arguments. In Sec.|3|we provide an overview of 
the general formalism for hidden variables network models 
and discuss how can we map the original SM into a model 
belonging this class of networks. In Sec. 0] we proceed to 
solve the mapped model, and provide analytical expres- 
sion for its main quantitative properties. The analytical 
results obtained are checked by means of direct numerical 



simulations of the orig inal SM in Sec. El Finally, in SecEl 
we draw the conclusions of our work. 



2 The static model 

The static model (SM) was introduced in Ref. as an 
algorithm to generate scale-free static (i.e. not growing) 
networks with any desired degree exponent 7 larger than 
or equal to 2. The model is defined as follows: We start 
from TV disconnected vertices, each one of them indexed 
by an integer number z, taking the values i — 1, . . . N. To 
each vertex, a normalized probability pi is assigned, given 
as function of the index i by 



1^.7 = 1 J 



(3) 



where a is a real number in the range a G [0,1]. The 
network is constructed iterating the following rules: Two 
different vertices i and j are randomly selected from the 
set of N vertices, with probability pi and pj, respectively. 
If there exists an edge between these two vertices, they are 
discarded and a new pair is randomly drawn. Otherwise, 
an edge is created between vertices i and j. This process is 
repeated until E — mN edges are created in the network, 
accounting for a fixed average degree (k) — 2E/N — 2m. 

This algorithm generates networks in which, by con- 
struction, there are no self-connections (a vertex joined to 
itself) not multiple connections (two vertices connected 
by more than one edge). The corresponding degree distri- 
bution can be estimated by means of a simple mean-field 
argument '25'. Since edges are connected to vertices with 
a probability given by the factor pi, we have that the prob- 
ability that any edge belongs to the vertex i, with degree 
fci, is given by 

In the large N limit, approximating sums by integrals, we 
have that, for < a < 1, 



N 



1 — a 



(5) 



Therefore, since J2j % = {k)N, we have from Eq. that 
h-PrJ2k, r^2m{l~a)f^) . (6) 



From this last expression, and using general arguments 
from network theory j^, we conclude that the degree 
distribution characterizing these networks has a scale-free 
form, P{k) k^'^ , with a degree exponent 



7=1 + -. 



(7) 



Thus, tuning the parameter a in the range [0, 1] it is pos- 
sible to generate networks with a degree exponent in the 
range 76 [2, 00]. 
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Just at this stage, it is possible to notice that the SM 
generates networks with built-in degree correlations. From 
Eq. © , we observe that the maximum degree, correspond- 
ing to the index i = 1, is given by 



h=i - 2to(1 - a)N". 



(8) 



This implies that the cut-ofF (or maximum expected de- 
gree) kc{N) in the network "5^^ scales with the network 
size as kc{N) ~ N". Now, it has been proved that, in or- 
der to have no correlations in the absence of multiple and 
self-connections, a scale-free networks with size N must 
have a cut-off scaling at most as ks{N) ^ N^^^ (the so- 
called structural cut-off) 3 1 . Therefore, the SM should 
yield correlated networks for values a > 1/2, i.e., for de- 
gree exponents in the interval 2 < 7 < 3, which correspond 
to those values empirically observed in real scale-free net- 
works. In the following Sections we will provide an ana- 
lytical description of the origin and form of these degree 
correlations. 



3 Mapping to a hidden variables network 
model 

In order to solve analytically the SM, it is useful to map 
it to a hidden variables network model |2()ll23ll24j . Hid- 
den variables network models are a generalization of the 
Erdos-Renyi model |32j in which vertices are assigned a 
tag (or hidden variable) whose statistical properties com- 
pletely determine the topological structure of the ensuing 
networks. 



3.1 General network models with hidden variables 

The class of network models with hidden variables is de- 
fined as follows 1201: Starting from a set of N disconnected 
vertices and a general hidden variable h, that can be a nat- 
ural or real number, we construct an undirected network 
with no self nor multiple connections, by applying these 
two rules: 

1. To each vertex i, a variable hi is assigned, drawn at 
random from the probability distribution p{h). 

2. For each pair of vertices i and j, with hidden variables 
hi and hj, respectively, an edge is created with prob- 
ability r{hi,hj) (the connection probability), where 
r{h,h') > is a symmetric function of h and h' . 

In this class of models, the degree distribution is given 



by 



(9) 



where the propagator g{k\h) gives the conditional proba- 
bility that a vertex with hidden variable h ends up con- 
nected to k vertices. The propagator is a normalized func- 
tion, '^f^ gik\h) = 1, whose generating function g(z\h), 
defined by 

g{z\h)^y^z'g{k\h), (10) 



fulfills in the general case the expression (201 
In h{z\h)^NY^ p{h') In [1 - (1 - z)r{h, h')] 



(11) 



Given the probabilities p{h) and r{h, h'), Eq. Hll|l must be 
solved and inverted in order to obtain the corresponding 
propagator and the degree distribution. Without solving 
this equation, however, we can still obtain some informa- 
tion on the connectivity properties of the network. Notic- 
ing that the first moment of g{k\h) is given by the first 
derivative of g{z\h), evaluated at z = 1, we that the av- 
erage degree of the vertices of hidden variable h, k{h), is 
given by 

k{h) ^Y.^g{k\h) = NY,p{h')r{h,h'), (12) 

k h' 

while the average degree takes the form 



(fc) = ^P(fc) = ^p(;j)fc(/j). 



(13) 



In order to characterize degree correlations in a gen- 
eral model with hidden variables, we need to provide an 
expression for the average degree of the neighbors of the 
vertices of degree k, knn{k)- Consider first the average de- 
gree of the neighbors of the vertices of hidden variable /i, 
knn{h). This quantity can be expressed as 



knnih) ^Y.~^{h')p{h'\h), 



(14) 



where p(h'\h) is the conditional probability that a ver- 
tex of hidden variable h is connected to a vertex of hidden 
variable h' . To compute this last quantity, we observe that 
the probability of drawing an edge from h to h' is propor- 
tional to the probability of finding an h' vertex, times the 
probability of creating an actual edge. Therefore, 



p{h'\h) 



p{h')r{h,h') 
Y.^„ p{h")r{h,h" 



Np{h')r{h,h') 
kih) ■ 



Thus, we have that 



knn(h) 



N 



Y,p{h')k(h')r{h,h'). 



(15) 



(16) 



Finally, the correlation function knn{k) can be shown to 
be given by [201 

fc„„(fc) = 1 + P{h)g{k\h)k^^{h). (17) 



3.2 Mapping the static model 

In order to map the SM into a hidden variables network 
model we need to provide a proper definition of the hid- 
den variables /i, their probability distribution p{h), and 
the connection probability r[h,h'). A natural choice for 
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the hidden variable is the index i associated to each ver- 
tex. On its turn, the connection probabihty r{i,j) can be 
defined as the probability that vertices i and j end up 
connected in the final network. With the original defini- 
tion of the SM, it is difficult to estimate this connection 
probability. In order to overcome this difficulty, we will 
consider a small variation of the algorithm defining the 
model. Within the original definition, in a first step of 
the model, a potential edge («, j) is selected, by randomly 
choosing a pair of vertices i and j, with probabilities pi 
and pj, respectively, as given by Eq. ©. In a second step, 
the potential edge is actually created if it did not 

exist previously, and this process is repeated until a given 
number of actual edges E = mN is reached, leading to a 
constant average degree (fc) = 2m. Thus, we can consider 
this as a microcanonical model, since the average degree 
is held fixed. This fact is in opposition with the spirit of 
hidden variables network models, in which the average de- 
gree is not constant, but tends to an asymptotic value for 
large network sizes [201 • We can place the SM within this 
network class by converting it to a canonical model, in 
which a fixed number E = mN of potential edges is cho- 
sen, and afterwards checked for their actual addition to 
the network. This canonical version of the SM will lead to 
a network with a number of edges smaller than or equal 
to E, and therefore to an average degree (k) < 2m. How- 
ever, we expect that this canonical version of the SM will 
coincide with the microcanonical original SM in the infi- 
nite network size limit, and to observe (k) 2m in the 
limit N (X. The good agreement between theoretical 
predictions derived from the first and simulations of the 
second will confirm this claim. 

Let us look at the edge creation process in the canon- 
ical version of the SM. If we allow for the possibility to 
choose a potential edge with i = j (self-connection), the 
probability of selecting the potential edge (i, j) is 2piPj 
if i j, and if i = j. In a more compact form, the 
probability of choosing the potential edge («, j) is 

7i"(i,j) = (2 - %)piPj, (18) 

where Sij is the Kronecker symbol. This probability is nat- 
urally normalized: If we sum 7r(i, j) over all the N{N+l)/2 
possible potential edges (including self-connections), we 
have 

i<j i<j i 

i<j i \ i / 

since the original distribution pi is normalized. 

The probability that, in the final network, the ver- 
tices i and j are connected is equal to the probability that 
the potential edge (z,j) has been selected at least once, 
which is the complementary probability that it has not 
been selected in the E trials made to generate the net- 
work. Therefore, for the canonical version of the SM, we 
have that the probability that vertices i and j are con- 



nected in the network is 

p,{t,j) = !-[!- 7r{t, j)f . (20) 

This expression can be further simplified by taking the 
limit of large N. We have that E — mN and pi ~ i~"{l — 
a)N"~^. Therefore, we can write 

Pc{^,J) = 1 - [1 - (2 - %)(1 - afN^'^-^z-^rT'' • 

.(21) 

Thus, in the limit N oo, we can approximate this ex- 
pression by an exponential, that yields the final result 

Pc{i,j) = 1 - exp [-(2 - %)to(1 - a)27V2"-ii-j-"] . 

(22) 

This is the probability that two vertices end up connected 
in the final network in the canonical version of the SM. 
Therefore, in the hidden variables version of the model we 
can set the connection probability 

r{i, j) = 1 - exp [-2to(1 - a)^N^°'^^i-°'j-°'] , (23) 

where we have neglected the Kronecker symbol, since in 
hidden variable models we do no allow for the possibil- 
ity of self-connections. A first conclusion can be extracted 
from this connections probability: it does not factorize in 
two independent functions of i and j. Therefore, degree 
correlations will be present in the model 20 . 

To complete the mapping, we finally need to give a pre- 
scription for the probability p{i) of a vertex having hidden 
variable (index) i. In the original definition of the model, 
the index is assigned deterministically to each vertex. Here 
we will assume an approximation already made for other 
models [30], that consists in considering the hidden vari- 
able i randomly assigned from the set {1,2,..., N}, with 
probability p{i) = 1/N. As we will see in the next Sec- 
tions, this assumption does not have a strong influence in 
most of the analytic results, when compared with numer- 
ical simulations of the original SM. 

4 Analytic solution 

4.1 Average degree 

Let us consider in the first place the behavior of the overall 
average degree, and the average degree of the vertices with 
index i. From Eq. H12|l . togheter with the definition of the 
probabilities p{i) and r{i,j), we have that 

N 

k{i)^NY,p{Mhj) 

N 

= ^ {1 - exp [-2m(l - af N^°'-\~°' J-"'] } . (24) 

Approximating sums by integrals, and performing the change 
of variables j = Nx, we are led to the expression 

k{i) ^N [ dx{l-exp [-2m(l - afN°'-^i-"x-°'] } . 

(25) 
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Since a < 1, the argument of the exponential is a de- 
creasing function of N. Therefore, in the hmit N —^ oo, 
we can perform a Taylor expansion of the integrand, and 
approximate 

k{i) f dx [2m(l - a)^N°'-H-"x-°'] (26) 



N- 



= 2m(l 



N 



(1 _iV"-i). 



(27) 



For large N, the last term in this expression tends to 1, 
and we recover the mean-field result obtained previously 
for the SM, Eq. ©. 

As for the average degree, we have from Eq. 



N 



(fc}=^p(^)fc(^) 



TV 



2to(1 - a) 



N 



N 



(1 -7V"-i) 



2m(l -iV"-^)^ 



(28) 

where again we have approximated sums by integrals. We 
observe that, for any finite network size, (k) < 2m. How- 
ever, in the thermodynamic limit N — > oo, we recover the 
fixed degree exponent (k) = 2m imposed by the SM. 



where we have approximated i as a continuous variable, 
performed the change of variables i = Nx~^^°', and ex- 
pressed kl — r{k + l), where r{z) is the standard Gamma 
function. The only dependence of expression on the net- 
work size is through the upper limit in the integral. There- 
fore, in the thermodynamic limit we can use the result |33| 



/oo 
e-''yy^dy = B-'-^r{l + A,B), 



(34) 



where r{z, a) is the incomplete Gamma function, to ob- 
tain 

_ [2m(l-a)]i/"r(fc-l/a,2m[l-a]) 

^^"^ - a r(fc + i) ■ ^^^^ 

In order to obtain the asymptotic behavior of the de- 
gree distribution for large fc, we note that r(z, a) r{z) 
for z — > oo. Therefore, for large k 

P(.).Mi^i^Ii^^^.--V.. (36) 
a r{k + 1) 

That is, we recover a scale- free degree distribution with 
a degree exponent 7 = l-|-l/a, as derived by mean- field 
arguments for the original SM. 



4.2 Degree distribution 

In order to compute the degree distribution, we must first 
solve Eq. for the generating function of the propaga- 
tor, g{z\i). For the probabilities p(i) and r{i, j) we are con- 
sidering, approximating sums by integrals and performing 
again the change of variables j = Nx, we have that 



hig{z\i) = N 



[ da;ln[l-(l- 



z X 



(29) 



X (l - exp {-2m(l - af N"-'^r°'x-"})] . (30) 

For a < 1, the argument in the exponential is again de- 
creasing in the large N limit. Therefore, expanding to first 
order the exponential, and then the logarithm inside the 
integral, we are led to 

In g{z\i) I dx{l - z)2m{l - af N°'-\-"'x-"' 

= {l-z)k{z). (31) 

Given Eq. (|31|l . we find that the propagator is finally given 
by a Poisson form: 

exp[-fc(i)]A:(i)'^ 



kl 



Knowing the form of the propagator, we can derive the 
degree distribution applying Eq. ©, i.e. 

i—\ 

X J dxexp[-2m{l- a)x]x''~^-'^^", (33) 



4.3 Degree correlations 

Next, we aim to calculate the average nearest neighbor 
degree of the vertices with degree k, knnik), in order to 
evaluate correlations. To do so, we first compute the av- 
erage nearest neighbor degree of the vertices with index 
i, knn{i), that is given by Eq. 112|l . Using the expression 
for k{i) that we have evaluated in Eq. H27|l in the large N 
limit, we have 

N 

~knn{i) = 5] [1 - exp {-2to(1 - a)2iV2a-l^-a^.-a|] 

i=i 

(37) 

We can proceed as usual, replacing sums by integrals. In 
this case, however, it is not possible to Taylor expand the 
integral after an appropriate change of variables, since the 
extra factor in the integral causes it to diverge in its 
lower limit. We must therefore keep the full exponential 
form. After some formal manipulations, we can write 



knn (*) 



i»{N^-» - 1) 



1-a 



(32) +_ / dxx-^^°'exp{-2m{l-afN^''~H-°'x} 



a ./I 



dxx-^l"" exp {-2to(1 - a)2iV"-ir"x} . 



After applying the identity Eq. (|34|l . we are led to the 
solution 



knn (0 



1) 



1-a 
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+ — [2m(l - a)27V2""ii-"]-i+i/" 



X {r(l- -, 2m(l - afN^"'-^i-°') 
a 



r{\ , 2m(l - afN°'~^i-"') 

a 



(38) 



As we will see in the following Section, the approxima- 
tion given by Eq. (|38|l is in fact not very good, and a much 
better agreement with numerical simulations is obtained 
by performing numerically the summation in the original 
discrete expression Eq. (|37|l . This fact is due to the effects 
of the continuum approximation in the index i, which are 
negligible at the level of the degree distribution, but show 
up at the level of correlations. 

Finally, for the average degree of the nearest neigh- 
bors of the vertices of degree fc, fc„„(fc), we resort to the 
expression Eq. (|17|l . taking the form for the SM 



knn{k) — 1 + 



1 



N 



NP{k)k 



^exp[-fc(^)]fc(^)'=fc„„(^). (39) 



This expression is far too complex to obtain even an asymp- 
totic expression in the continuous i approximation, so we 
will compare numerical simulations with a direct numeri- 
cal evaluation of the summation in Eq. (|39|l . 



5 Numerical simulations 

We have checked the analytical predictions presented in 
the previous Section by means of extensive numerical sim- 
ulations of the original SM. We have generated networks 
with a variable, m = 3 and size N = 10^. All results are 
averaged over 10"^ realizations for each value of the pa- 
rameter a. Simulation were performed as follows: At each 
iteration, we extract a pair of real numbers according to a 
power-law probability distribution with exponent a, nor- 
malized between 0.5 and N + 0.5. Number are extracted 
using the Monte Carlo inversion method [331 ■ Then, we 
approximate each number to the nearest integer, so that 
the resulting pair is composed by integers between 1 and 
N. These are the two candidate vertices to be connected 
by and edge. If the proposed pair is composed by two 
identical numbers, or it has been extracted before, the ex- 
traction is rejected and repeated until two valid vertices 
are proposed. We iterate this procedure until a network 
oi E = mN edges is created. This algorithm corresponds 
exactly to the original SM. The only modification is that 
the probability distribution according to which we extract 
the candidate edges is not discrete, but continuous. Any- 
way, it is possible to see that the results of the proposed 
procedure are indistinguishable from those obtained from 
methods that start directly from a discretized distribu- 
tion, but require more computation time (for example, by 
using the rejection method [33j). 

In Fig. n we plot the average degree of the vertices with 
index i, k(i) for two different values of a, namely a ~ 0.55 
and a = 0.8, which correspond to the degree exponents 




Fig. 1. Average degree of the vertices with index i in the orig- 
inal SM, for two different values of a. The solid lines represent 
the theoretical value given by Eq. l|?7|l . 




Fig. 2. Degree distribution in the SM for two different values of 
a. The dotted and dashed lines represent the theoretical value 
given by Eq. . In the inset: enlargement of the tail of the 
distribution. Peaks due to the discretization of the degree are 
visible. Dashed vertical lines represent the theoretical values of 
the centers of such peaks. 



7 — 2.82 and 7 = 2.25, respectively. In both cases, the an- 
alytical result, as as given by Eq. 1)27(1 . fit almost perfectly 
the curves emerging from numerical simulation. The same 
happens for the degree distribution, shown in Fig. |2 for 
the two values of a considered. As we can see from this 
Figure, the complete expression calculated in Eq. H35|l fits 
exactly the whole distribution, except at very large val- 
ues of k. This discrepancy, due to the finite size of the 
networks, is easy to understand. From Ea. (|27() . we can 
observe that large values of k correspond to small values 
of the index i. In this region, the continuous i and k ap- 
proximation made in all calculations is expected to fail, 
and the index i to show its true discrete nature. Indeed, 
this fact can be clearly observed in the inset in Fig. (21 
where we plot a close-up of the tail of the degree distribu- 
tion obtained for a — 0.8, obtained from averaging over 
10'^ network samples. This plot shows a set of peaks, cor- 
responding to the first values of the index i, from 1 to 
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Fig. 3. Average nearest neighbor degree of the vertices with Fig. 4. Average nearest neighbor degree of the vertices with 



index i in the SM for two different values of a. The solid lines 
represent the theoretical value given by the numerical summa- 
tion of Eq. 1371 . The dashed lines correspond to the analytical 
approximation in the Eq. I|88^ . 



degree k, in the SM for two different values of a. The dashed 
lines correspond to the numerical summation of Eq. 1391 . 
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8. The centers of the peaks are well approximated by the 
analytical k{i) function given in Eq. (|27|l . and represented 
by means of vertical dotted hnes, except for very small 
values of i. The width of the peaks is accounted for by 
the fluctuations in the value of k in the different network 
samples. 

In Fig. O we report the average nearest neighbors de- 
gree of the vertices with index i. The dashed line repre- 
sents the theoretical approximation obtained in Eq. (|38l) . 
We find a percentually small difference between calcula- 
tion and simulation. This difference can be attributed to 
the effect of the continuous approximation. Indeed, if we 
numerically calculate the sum of Eq. l|57|l and report it 
in the plot (continuous line), we obtain a better fit of 
the simulation results. A good agreement between theory 
and simulation is obtained as well in the plot of the av- 
erage nearest neighbor degree of the vertices with degree 
k, Fig. 01 at least for sufficiently large values of a. Here 
theoretical value is obtained directly from numerical sum- 
mation of Eq. H39|l , that cannot be approximated analyti- 
cally in a simple way. The correlation function displays an 
almost constant behavior at low degrees and a decreasing 
slope at high degrees, i.e. a regime without any correla- 
tion followed by one characterized by strong disassortative 
mixing. The emergence of these correlation is connected 
to the absence of multiple and self-connections. By doing 
this, we bias the natural tendency of high degree vertices 
to have some connections into each other, favoring their 
linking to small degree vertices, and, therefore, generat- 
ing negative correlations in the degree. This phenomenon 
appears to be extremely relevant for small values of 7. 
In the case a = 0.8, for example, we can observe in the 
average neighbor connectivity a decay of more than one 
decade in about two decades of the degree. On the other 
hand, the analytical solutions does not behave so well for 
large values of 7 (small a), probably due to the accumu- 
lated effect of all the approximations made in obtaining 
this expression. 



In this work we have presented an analytic solution of the 
static model which has been recently proposed as a 
paradigmatic scale-free non-growing network model. The 
solution is obtained via a mapping of the SM into a hid- 
den variables network model, that represents its canonical 
counterpart, i.e. in which the number of edges is not held 
fixed, but whose average degree tends to a constant in the 
infinite network size limit. We have derived analytically 
the properties of the mapped hidden variables network 
model and checked the predictions by means of exten- 
sive numerical simulations of the original model. The good 
agreement observed implies that the canonical version of 
the SM is identical to the original version in the ther- 
modynamic limit — !■ 00. It is particularly noteworthy 
that our analytical calculations have allowed us to eval- 
uate the correlations induced in the model by the physi- 
cal condition of absence of self and multiple connections. 
The detected amount of correlation is considerable and 
can have a strong influence both on the topology of the 
networks, and on dynamics running on top of them. The 
presence of these correlations thus casts some shadows on 
the usefulness of the SM as a benchmark for mean-field so- 
lutions of dynamical processes, which are usually obtained 
in the uncorrelated limit. In this sense, the recently pro- 
posed Uncorrelated Configuration Model |35| appears to 
be a more adequate instrument for the investigation of the 
effects of scale-invariance in network topology and dynam- 
ics, without the perturbations induced by the presence of 
correlations. 
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